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Abstract 

In the alignment limit of a multi-doublet Higgs sector, one of the Higgs mass eigenstates 
aligns with the direction of the scalar field vacuum expectation values, and its couplings 
approach those of the Standard Model (SM) Higgs boson. We consider CP-conserving 
Two-Higgs-Doublet Models (2HDMs) of Type I and Type H near the alignment limit in 
which the lighter of the two CP-even Higgs bosons, h, is the SM-like state observed at 
125 GeV. In particular, we focus on the 2HDM parameter regime where the coupling 
of h to gauge bosons approaches that of the SM. We review the theoretical structure 
and analyze the phenomenological implications of the regime of alignment limit without 
decoupling, in which the other Higgs scalar masses are not significantly larger than 
and thus do not decouple from the effective theory at the electroweak scale. For the nu¬ 
merical analysis, we perform scans of the 2HDM parameter space employing the software 
packages 2HDMC and Lilith, taking into account all relevant pre-LHC constraints, the 
latest constraints from the measurements of the 125 GeV Higgs signal at the LHC, as well 
as the most recent limits coming from searches for heavy Higgs-like states. We contrast 
these results with the alignment limit achieved via the decoupling of heavier scalar states, 
where h is the only light Higgs scalar. Implications for Run 2 at the LHC, including 
expectations for observing the other scalar states, are also discussed. 
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1 Introduction 


The minimal version of the Standard Model (SM) contains one complex Higgs doublet, 
resulting in one physical neutral CP-even Higgs boson after electroweak symmetry breaking. 
The discovery [1,2] of a new particle with mass of about 125 GeV [3] and properties that match 
very well those expected for a SM Higgs boson was a real triumph of Run 1 of the LHC. Fits of 
the Higgs couplings performed by ATLAS [4] and CMS [5] show no significant deviations from 
SM expectations. (A combined global ht of the Higgs couplings based on the Run 1 results was 
performed by some of us in [6].) However, one has to keep in mind that the present precisions 
on the Higgs couplings are, roughly, of the order of tens of percent, so substantial deviations are 
still possible. Indeed, the SM is not necessarily the ultimate theoretical structure responsible for 
electroweak symmetry breaking, and theories that go beyond the SM, such as supersymmetry, 
typically require an extended Higgs sector [7-10]. Hence, the challenge for Run 2 of the LHC, 
and other future collider programs, is to determine whether the observed state is the SM Higgs 
boson, or whether it is part of a non-minimal Higgs sector of a more fundamental theory. 

In this paper, we take Two-Higgs-Doublet Models (2HDMs) of Type I and Type H [11] 
as the prototypes for studying the effects of an extended Higgs sector. Our focus will be on 
a particularly interesting limit of these models, namely the case in which one of the neutral 
Higgs mass eigenstates is approximately aligned with the direction of the scalar held vacuum 
expectation values. In this case, the coupling to gauge bosons of the Higgs boson observed 
at the LHC tends towards the SM limit, Cy —> 1.^ This so-called alignment limit is most 
easily attained in the decoupling limit [12], where all the other non-SM-like Higgs scalars of 
the model are heavy. However, the alignment limit of the 2HDM can also be achieved in a 
parameter regime in which one or more of the non-SM-like Higgs scalars are light (and in some 
cases very light). This region of alignment without decoupling is a primary focus of this paper. 

An extensive review of the status of 2HDMs of Type I and Type H was given in [13,14]. 
Interpretations of the recently discovered Higgs boson at 125 GeV in the context of the 2HDMs 
were also studied in [15-21]. The possibility of alignment without decoupling was first noted in 
[12] and further clarified in [22,23]. Previous studies of alignment without decoupling scenarios 
in the light of the LHG Higgs results were conducted in [24-26] . The specific case of additional 
light Higgs states in 2HDMs with mass below ~ 125/2 GeV was studied in [27]. 

Considering experimental as well as theoretical uncertainties, the expected precision for 
coupling measurements at the LHC after collecting 300 fb“^ of data is about 4-6% for the 
coupling to gauge bosons, and of the level of 6-13% for the couplings to fermions [28]. The 
precision improves by roughly a factor of 2 for at the high-luminosity run of the LHC with 
3000 fb“L At a future e+e" international linear collider (ILC) with y/s = 250 GeV to 1 TeV, 
one may measure the couplings to fermions at the percent level, and the coupling to gauge 
bosons at the sub-percent level. A detailed discussion of the prospects of various future colliders 
can be found in [28]. 

We take this envisaged ~ 1% accuracy on Cy as the starting point for the numerical analysis 
of the alignment case. Concretely, we investigate the parameter spaces of the 2HDMs of Type I 
and Type H assuming that the observed 125 GeV state is the h, the lighter of the two GP-even 

^We use the notation of coupling scale factors, or reduced couplings, employed in [6]; Cy {V = W,Z) for 
the coupling to gauge bosons, Cu,d for the couplings to up-type and down-type fermions and C-y^g for the 
loop-induced couplings to photons and gluons. 
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Higgs bosons in these models, and imposing that Cy > 0.99 (note that \Cv\ < 1 in any model 
whose Higgs sector consists of only doublets and/or singlets). The case of the heavier CP-even 
H being the state at 125 GeV is discussed in a separate paper [29]. 

Taking into account all relevant theoretical and phenomenological constraints, including the 
signal strengths of the observed Higgs boson, as well as the most recent limits from the non¬ 
observation of any other Higgs-like states, we then analyse the phenomenological consequences 
of this scenario. In particular, we study the variations in the couplings to fermions and in the 
triple-Higgs couplings that are possible as a function of the amount of alignment when the other 
Higgs states are light, and contrast this to what happens in the decoupling regime. Moreover, 
we study the prospects to discover the additional Higgs states when they are light. 

The public tools used in this study include 2HDMC [30] for computing couplings and decay 
widths and for testing theoretical constraints within the 2HDM context, Lilith 1.1.2 [31] for 
evaluating the Higgs signal strength constraints, and SusHi-1 .3.0 [32] and VBFNLO-2 .6.3 [33] 
for computing production cross sections at the LHC. 

The paper is organised as follows. In Section 2 we first review the theoretical structure of 
the 2HDM. A softly-broken discrete Z 2 -symmetric scalar potential is introduced using a basis of 
scalar doublet fields (called the Z 2 -basis) in which a the symmetry is manifest. The Higgs basis 
is then introduced, which provides an elegant framework for exhibiting the alignment limit. We 
then provide a comprehensive discussion of the Higgs couplings in the alignment regime. In 
Section 3, we explain the setup of the numerical analysis and the tools used. The results are 
presented in Section 4. Section 5 contains our conclusions. In Appendix A, detailed formulae 
relating the quartic coefficients of the Higgs potential in the Z 2 -basis to those of the Higgs 
basis are given. Some useful analytical expressions regarding the trilinear Higgs self-couplings 
in terms of physical Higgs masses are collected in Appendix B. 

NB: This version of the paper has been updated to include exactly the same 
constraints as Part 2 with mn = 125 GeV [29]. 

2 CP-conserving 2HDM of Types I and II 

In this section, we review the theoretical structure of the two-Higgs doublet model. Compre¬ 
hensive reviews of the model can also be found in, e.g., [12,23,34,35]. In order to avoid tree-level 
Higgs-mediated flavor changing neutral currents (FCNCs), we shall impose a Type-I or H struc¬ 
ture on the Higgs-fermion interactions. This structure can be naturally implemented [36,37] by 
imposing a discrete Z 2 symmetry on the dimension-four terms of the Higgs Lagrangian. This 
discrete symmetry is softly-broken by mass terms that appear in the Higgs scalar potential. 
Nevertheless, the absence of tree-level Higgs-mediated FCNCs is maintained, and FCNC effects 
generated at one loop are all small enough to be consistent with phenomenological constraints 
over a significant fraction of the 2HDM parameter space [38-41] . 

Even with the imposition of the softly-broken discrete Z 2 symmetry mentioned above, new 
CP-violating phenomena in the Higgs sector are still possible, either explicitly due to a physical 
complex phase that cannot be removed from the scalar potential parameters or spontaneously 
due to a CP-violating vacuum state. To simplify the analysis in this paper, we shall assume 
that these CP-violating effects are absent, in which case one can choose a basis of scalar doublet 
Higgs fields such that all scalar potential parameters and the two neutral Higgs field vacuum 
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expectation values are simultaneously real. Moreover, we assume that only the neutral Higgs 
helds acquire non-zero vacuum expectation values, i.e. the scalar potential does not admit the 
possibility of stable charge-breaking minima [43,44], 

We hrst exhibit the Higgs scalar potential, the corresponding Higgs scalar spectrum and the 
Higgs-fermion interactions subject to the restrictions discussed above. Motivated by the Higgs 
data, we then examine the conditions that yield an approximately SM-like Higgs boson. 


2.1 Higgs scalar potential 

Let $1 and $2 denote two complex Y = 1, 811(2)^;^ doublet scalar fields. The most general 
gauge invariant renormalizable scalar potential is given by 

V = - K 2 d>ld >2 + h.c.j + |Ai(d>ld>i )2 + + A 3 (l>Il>i)(d>^d> 2 ) 


+A4(l>ll>2)(d>^d>i) + 


iA5(d>ld>2)'+ [A6(l>Il>i) + A7(l>^l>2)]l>ld>2 + h.c.} . (1) 


In general, m^ 2 ) -^ 5 ) ^6 and A 7 can be complex. As noted above, to avoid tree-level Higgs- 
mediated FCNCs, we impose a softly-broken discrete Z 2 symmetry, $1 ^ -|-$i and $2 ^ —$2 
on the quartic terms of Eq. (1), which implies that Ae = A 7 = 0, whereas m ^2 0 is allowed. 

In this basis of scalar doublet fields (denoted as the Z 2 -basis), the discrete Z 2 symmetry of 
the quartic terms of Eq. ( 1 ) is manifest. Eurthermore, we assume that the scalar fields can be 
rephased such that m ^2 and A 5 are both real. The resulting scalar potential is then explicitly 
CP-conserving. 

The scalar fields will develop non-zero vacuum expectation values if the Higgs mass matrix 

has at least one negative eigenvalue. We assume that the parameters of the scalar potential 
are chosen such that the minimum of the scalar potential respects the U(l)gjy[ gauge symmetry. 
Then, the scalar field vacuum expectations values are of the form 


(^ 1 ) 


1 

72 


0 

Vl 


(^ 2 ) 


1 

72 



( 2 ) 


As noted in Appendix B of Ref. [ 12 ], if |m^ 2 l ^ -^ 5 bill^ 2 !, then the vacuum is CP-conserving 
and the vacuum expectation values vi and V 2 can be chosen to be non-negative without loss of 
generality. In this case, the corresponding potential minimum conditions are:^ 


"^11 = ”^12^/3 - 7“^ + ^3455^) , (3) 

^22 = ^12^7 ~ + A 345 CJ) , (4) 

where we have defined: 


-^345 = A3 + A4 + A5 , 


where 0 < /5 < and 


2 2,2 

V = VI + V2 


9^ 


V2 


tjs = tan /3 = — , 


Vl 


(246 GeVf . 


( 5 ) 

( 6 ) 


^Here and in the following, we use the shorthand notation C/3 = cos/ 3 , S/3 = sin/ 3 , = cos a, = sin a, 
C2/3 = cos 2/d, S2/3 = sin 2 / 3 , = cos (/3 — a), = sin(/d — a), etc. 
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Of the original eight scalar degrees of freedom, three Goldstone bosons [G^ and G) are 
absorbed (“eaten”) by the and Z. The remaining five physical Higgs particles are: two 
CP-even scalars {h and H, with < niH), one CP-odd scalar (^4) and a charged Higgs pair 
The resulting squared-masses for the CP-odd and charged Higgs states are 

m\ = , (7) 

m^± =m\ + |n^(A 5 - A 4 ), ( 8 ) 


where 


—2 _ 


m = 


'S 2/3 


( 9 ) 


The two neutral CP-even Higgs states mix according to the following squared-mass matrix: 

7 ' 77 ^ ' 777 ^ ^''77 . 77 ' '7 77 ^ 1 - (lo) 


AiU^c^ + {m\ + A 5 U^)s| [A 345 W^ - {rr7 + A 5 U^)]s/ 3 C /3 

[A 345 t'^ - {rn\ + Asu^)] X 2 v‘^s'j) + {m\ + 


Defining the physical mass eigenstates 

H = (v^Red*® — vi)ca + (v^Red*® — V2)sa , 
h = —(v^Red*® — ui)sq, -|- (\/2Red>2 — V2)ca , 

the masses and mixing angle a are found from the diagonalization process 


(11) 

( 12 ) 


m 


H 


0 

0 ml 


Srv Cn 


M7 Ml 


M\^M 


Cry Sr 


Mlicl + 2Ml2CaSa + MI 2 SI M^i^ “ (A4^2 - Mh)^aCo 

2^12(4 - 4 ) + (-A^22 - M 7 )^aCa “ ^M^CaSa + M %7 


(13) 


^I2'^cx^cx 


^22^0i 


ml 


Note that the two equations, Tr MT' — m^j + ml and det M^ = mj^ml, yield the following 
result: 

' ' ( 14 ) 

(15) 


1^7 = ijn]j - Ml7{Ml^ - ml) = ^J{Ml2-ml){M 
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Explicitly, the squared-masses of the neutral CP-even Higgs bosons are given by 

77,h ~ \ \/ 7 \l + M\2 =t a] , 

where < mn and the non-negative quantity A is defined by 




(16) 


The mixing angle ck, which is defined modulo tt, is evaluated by setting the off-diagonal 
elements of the CP-even scalar squared-mass matrix given in Eq. (13) to zero. It is often 
convenient to restrict the range of the mixing angle to |q;| < Itt. In this case, is non-negative 
and is given by 


Cry 


A + M?, -M^ 


22 


2 A 


lM7-ml 

9 9 ’ 

mjj-mi 


(17) 
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and the sign of is given by the sign of Aif 2 - Explicitly, we have 




V2MI2 

v/A(A + Ml^- MI2) 


sgn(Af?2) 




ml - Ml, 

ml - ml 


( 18 ) 


In deriving Eqs. (17) and (18), we have assumed that 7 ^ mn- The case of ruh = ran is 
singular; in this case, the angle a is undefined since any two linearly independent combinations 
of h and H can serve as the physical states. In the rest of this paper, we shall not consider this 
mass-degenerate case further. 


2.2 SM-limit in the Higgs basis 

The scalar potential given in Eq. (1) is expressed in the Z 2 -basis of scalar doublet helds in 
which the Z 2 discrete symmetry of the quartic terms is manifest. It will prove convenient to 
re-express the scalar doublet fields in the Higgs basis [45,46], defined by 

H, = , H 2 = (^^0 = + $2C(3 , (19) 

SO that (i7°) = vl\f2 and (77°) = 0- The scalar doublet i7i possesses SM tree-level couplings 
to all the SM particles. Therefore, if one of the CP-even neutral Higgs mass eigenstates is 
SM-like, then it must be approximately aligned with the real part of the neutral field i7°. 

The scalar potential, when expressed in terms of the doublet fields, and 772 , has the 
same form as Eq. (1), 

V = Ti77|77i + ^2771772 + T 3 [ 77 l 772 + h.c.j + + Z^{E\E^){H\H< 2 ) 

+Z 4 ( 77 | 772 )( 772 ^ 77 i) + |lZ 5 ( 77 | 772 )' + {Z^{E\Ex) + Z^{hIH 2 )\h\H 2 + h.c.} , ( 20 ) 
where the Yi are real linear combinations of the and the Zi are real linear combinations of 


the Aj. In particular, since Ae = A7 = 0 , we have [ 46 , 47 ]^ 

Zi = AiC^ -|- A2S^ -|- |A 345 S 2 ,g , (21) 

Z2 = Ais^ -|- A2C^ + 7A345S2/3, (22) 

= i^2/9 [-^1 + A2 — 2A345] -|- Aj, (for z = 3,4 or 5) , ( 23 ) 

Zq = “ ^2S^ — A345C2/3] , ( 24 ) 

Z^ = — |s 2/3 [AiS^ — A2C^ + A345C2/3] . ( 25 ) 

Since there are five nonzero A* and seven nonzero Zj, there must be two relations. The following 
two identities are satisfied if /3 7^ 0 , ^tt, Itt [ 47 ] d 

Z2 = Z\ Y 2(.^6 -|- Z-j) cot 2/3 , ( 26 ) 

.^345 = Zi Y ‘^Zq cot 2/3 — {^Zq — Zt) tan 2/3 , ( 27 ) 


^To make contact with the notation of Ref. [12], A = Zi, Ay = Z2, At = .^3 -h -^4 — -^5, At = — -^4, 

= Z\ — A = — Zq and Xjj = — Zj. 

^For /d = 0 , ^TT, the Z2-basis and the Higgs basis coincide, in which case Zq = Zj = ^ and Zi, Z2, Z345 are 
independent quantities. For (3 = the two relations are Zi = Z2 and Zq = Z^^ and Z345 is an independent 
quantity. 
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where Z 345 = ^3 + ^4 + Z 5 . One can invert the expressions given in Eqs. (21)-(25), subject to 
the relations given by Eqs. (26) and (27). The results are given in Appendix A. 


The squared mass parameters are given by 

Yi = + ml^sl - , (28) 

¥2 = + m22C^ + mi2S2/3, (29) 

^3 = - mii)s2/3 - mi2C2/3. (30) 

Yi and Y 3 are fixed by the scalar potential minimum conditions, 

Yi = -lZiv\ Ys = -lZev\ (31) 

Using Eqs. (9) and (31), we can express fn^ in terms of 12, Zi and Zq, 

= >2 + + Zqv'^ cot 2(3 . (32) 

The masses of and A are given by 

m^± = y2 + \Z^v^ , (33) 

= I 2 + I(-^3 Z 4 ^ — Z^'iP '. (34) 


It is straightforward to compute the CP-even Higgs squared-mass matrix in the Higgs ba¬ 
sis [45,48], 


M 


2 _ 

H — 


(Ziv^ 

[Zev^ 


Zev^ \ 
m\ + Z^v^J 


(35) 


From Eq. (35), one can immediately derive the conditions that yield a SM-like Higgs boson. 
Since {H^) = n/\/2 and (Tf®) = 0; couplings of Hi are precisely those of the Standard 
Model. Thus a SM-like Higgs boson exists if \/2 Re — n is an approximate mass eigenstate. 
That is, the mixing of and 77° subdominant, which implies that either \Zq\ <C 1 and/or 
m\-\-Z^v^ ^ Ziv'^, Zqv"^. Moreover, if in addition Ziv^ < m^ + Z^v"^, then h is SM-like, whereas 
if Ziv'^ > m\ + Z^v'^, then H is SM-like. In both cases, the squared-mass of the SM-like Higgs 
boson is approximately equal to Ziv'^. 

The physical mass eigenstates are identified from Eq. (11), (12) and (19) as 



H = (v^ReTf® — v)ci3-a — V2Re , 

(36) 


h = {'/2ReHi — v)si3-a + v^Re ■ 

(37) 

Then, Eqs. 

(15) and (16) yield 



^H,h = l['n^A + {Zi + Z5)v‘^ ± Ah] , 

(38) 

where 

Ah^^ K + (^5 - Zi)v ^]' + 4Z|n4 . 

(39) 

In addition. 

, Eq. (14) yields 



{m% - Ziv'^){Ziv'^ - ml) . 

(40) 
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Comparing Eqs. (11) and (12) with Eqs. (36) and (37), we identify the corresponding mixing 
angle by a — /3, which is defined modulo tt. Diagonalizing the squared mass matrix, Eq. (35), 
it is straightforward to derive the following expressions: 


ry ‘2 2 2 1 2 2 

ZlV = , 

(41) 

ZqV^ = {ml - m%)s^_aC/3-a , 

(42) 

2 1 V 2 2 2 1 2 2 

m^ + Z^v = . 

(43) 

It follows that 


ml= (zi + Zq , 

V ^/3-a J 

(44) 

rriH = + Zq ]v. 

V S(3-aJ 

(45) 

Note that Eq. (42) implies that^ 


^6^(3—a^(3—a — 0 . 

(46) 


One can also derive expressions for and either directly from Eqs. (41) and (42) or 
by using Eqs. (17) and (18) with a replaced by a — /3. Using Eq. (46), the sign of the product 
Sis-aCp-a is fixed by the sign of Z^. However, since /? — a is defined only modulo tt, we are 
free to choose a convention where either or is always non-negative.® In a convention 
where s^-a is non-negative (this is a convenient choice when the h is SM-like), 


^(3—a 


sgn(Z6) 


/ 


Ziv^^ — ml 
m%-ml 


_ -^ 6 ^^ _ 

- Ziv"^) ’ 


(47) 


where we have used Eq. (40) to obtain the second form for in Eq. (47). 

Finally, we record the following useful formula that is easily obtained from Eqs. (7) and 
(A.IO),^ 

= m\ -\- Z^v^ + ^(Zq — Z 7 )v^ tan 2/3 . (48) 

Combining Eq. (48) with Eqs. (42) and (43) yields 

Zrv^ = {ml - m]j)sp-ccCp-a + 2 cot 2/3 + mlcl_^ - m^] . (49) 

Using Eqs. (26) and (27), one can likewise obtain expressions for ^ 2 ^^ and ^ 345 ^;^ in terms of 
and However, these expressions are not particularly illuminating, so we do not 
write them out explicitly here. 

^Having established a convention where 0 < < ^tt, we are no longer free to redehne the Higgs basis held 

H 2 —H 2 . Consequently, the sign of Zq is meaningful in this convention. 

^Such a convention, if adopted, would replace the convention employed in Eq. (17) in which is taken to 
be non-negative. 

^In Eq. (48), the term in the expression for m ^ that is proportional to {Zq — Zj)v‘^ tan 2/3 is never greater than 
0{v‘^) for all values of P, since Eqs. (24) and (25) imply that {Zq — Z^) tan2^5 = — A 2 — 2 A 345 ) < 0(1). 
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2.3 Higgs couplings and the alignment limit 

As noted in the previons subsection, the Higgs basis field Hi behaves precisely as the 
Standard Model Higgs boson. Thus, if one of the neutral CP-even Higgs mass eigenstates is 
approximately aligned with \/2Rei/° — v, then its properties will approximately coincide with 
those of the SM Higgs boson. Thus, we shall define the alignment limit as the limit in which the 
one of the two neutral CP-even Higgs mass eigenstates aligns with the direction of the scalar 
field vacuum expectation values. Defined in this way, it is clear that the alignment limit is 
independent of the choice of basis for the two Higgs doublet fields. Nevertheless, the alignment 
limit is most clearly exhibited in the Higgs basis. In light of Eqs. (36) and (37), the alignment 
limit corresponds either to the limit of ^ 0 if h is identified as the SM-like Higgs boson, 
or to the limit of sg^a 0 ii H is identified as the SM-like Higgs boson. 

Consider first the case of a SM-like h, with ~ 125 GeV. In this case, Ziv'^ < m\ + Z^v'^, 
\ci 3 -a\ -C 1, and ml ~ Ziv"^. It follows from Eq. (47) that the alignment limit can be achieved 
in two ways: (i) Ze —?• 0 or (ii) mn v. The case of mn ^ v (or equivalently Y 2 ^ v) is 
called the decoupling limit in the literature.® In this case, one finds that mn ~ toa ^ ttih^ 1 so 
one can integrate out the heavy scalar states below the scale of mu- The effective Higgs theory 
below the scale mn is a theory with one Higgs doublet and corresponds to the Higgs sector of 
the Standard Model. Thus not surprisingly, /j is a SM-like Higgs boson. However, it is possible 
to achieve the alignment limit even if the masses of all scalar states are similar in magnitude 
in the limit of Zq 0. This is the case of alignment without decoupling and the main focus of 
this study. Einally, if both \Zq\ <C 1 and mn TUh are satisfied, the alignment is even more 
pronounced; when relevant we shall denote this case as the double decoupling limit. 

For completeness we note that in the case of a SM-like H we have Ziv"^ > m\ + Z^v'^, 
|s/ 3 _q,| -C 1 and ~ Ziv'^. Here, it is more convenient to employ a convention where cg-a is 
non-negative. One can then use Eqs. (40), (46) and (47) to obtain an expression for s^_a. In 
a convention where cg^a is non-negative. 


^(3—a 


sgn(Z6) 


/ 


m%-ml 


_ -Zev^ _ 

y/{mjj-ml){Ziv‘^-ml) 


(50) 


Taking uih ~ 125 GeV, there is no decoupling limit as in the case of a SM-like h. However, 
the alignment limit without decoupling can be achieved in the limit oi Zq ^ 0. This case will 
be discussed in detail in [29]. 

We now turn to the tree-level Higgs couplings. Denoting the SM Higgs boson by hsu^ the 
coupling of the CP-even Higgs bosons to VV (where V = or Z) normalized to the 
coupling is given by 

Cy = • (51) 

As expected, if /i is a SM-like Higgs boson then Cy 1 in the alignment limit, whereas if H is 
a SM-like Higgs boson then Cy 1 in the alignment limit. 

^More precisely, we are assuming that ^ \Zq\v‘^. Since Zq is a dimensionless coefficient in the Higgs 
basis scalar potential, we are implicitly assuming that Zq cannot get too large without spoiling perturbativity 
and/or unitarity. One might roughly expect \Zq\ < 47r, in which case mn ^ v provides a reasonable indication 
of the domain of the decoupling limit. 







Table 1: Four possible Z 2 charge assignments that forbid tree-level Higgs-mediated FCNC 
effects in the 2HDM [49]. 



d>i 

$2 



^R 

tL, bL, I'Ll Cl 

Type I 

+ 

— 

— 

— 

— 

+ 

Type H 

+ 

— 

— 

+ 

+ 

+ 

Type X (lepton speciflc) 

+ 

— 

— 

— 

+ 

+ 

Type Y (flipped) 

+ 

— 

— 

+ 

— 

+ 


Next, we consider the Higgs boson couplings to fermions. The most general renormalizable 
Yukawa couplings of the two Higgs doublets to a single generation of up and down-type quarks 
and leptons (using third generation notation) is given by 

-Fh.c., 
(52) 

where ei 2 = —621 = 1 , en = 622 = 0 , Ql = {ti ) and Ll = ) ^l) are the doublet 

left handed quark and lepton fields and Ir, bn and cr are the singlet right-handed quark and 
lepton fields. However, if all terms in Eq. (52) are present, then tree-level Higgs-mediated 
FCNCs would be present, in conflict with experimental constraints. To avoid tree-level Higgs- 
mediated FCNCs, we extend the discrete Z 2 symmetry to the Higgs-fermion Lagrangian. There 
are four possible choices for the transformation properties of the fermions with respect to Z 2 , 
which we exhibit in Table 1. 

For simplicity, we shall assume in this paper that the pattern of the Higgs couplings to down- 
type quarks and leptons is the same. This leaves two possible choices for the Higgs-fermion 
couplings [ 11 ]: 


-^Yuk = [yitRQ'M+yrtRQiH] 


Type I: = 0 , (53) 

Typell: y^ = y^ = y^^ = 0 . (54) 


In particular, the pattern of fermion couplings to the neutral Higgs bosons in the Type I and 
Type H models is exhibited in Table 2 . 

In the strict alignment limit, the fermion couplings to the SM-like Higgs boson should 
approach their Standard Model values. To see this explicitly, we note the identities. 


cos a 


sin 

/5 

sin 

a 

cos 

/5 

sin 

a 

sin 


cos 

a 


cos 13 


^iS-a + cot ) 

Sjj—oc tan (3cij—a , 
Cfj—d, cot f3si^—a , 
—cx T tan [3sp —q. . 


(55) 

(56) 

(57) 

(58) 


If h is the SM-like Higgs boson, then in the limit of Cjs^a ^ 0, the fermion couplings of h ap¬ 
proach their Standard Model values. However, if tan /5 ^ 1 , then the alignment limit is realized 
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Table 2 : Tree-level vector boson couplings Cy {V = W,Z) and fermionic couplings Cp normal¬ 
ized to their SM values for the two scalars h, H and the pseudoscalar A in Type I and Type II 
Two-Higgs-doublet models. 



Type I and H 

Type I 

Type H 

Higgs 

VV 

up quarks 

down quarks 
and leptons 

up quarks 

down quarks 
and leptons 

h 

sin(/5 — a) 

coscr/ sin/3 

cos a/ sin (3 

cos a/ sin/3 

—sin a/ cos /3 

H 

cos{P — a) 

sina/ sin/3 

sin a/ sin (3 

sin a/ sin/3 

cos a/ cos /3 

A 

0 

cot (3 

— cot (3 

cot (3 

tan /3 


in the Type-II Yukawa couplings to down-type fermions only if |c; 3 _a|tan /3 <C 1. That is, if 
\ci 3 -a\ ^ 1 but \ci 3 -a \ tan/5 ~ 0{1), then the hVV couplings and the hit couplings are SM-like 
whereas the hbb and hT^T~ couplings deviate from their Standard Model values. Thus the 
approach to the alignment limit is delayed when tan ^ 1. We denote this phenomenon as the 
delayed alignment limit. Similar considerations apply if cot/5 ^ 1 ; however, this region of pa¬ 
rameter space is disfavored as the corresponding hti coupling quickly becomes non-perturbative 
if cot 13 is too large. 

Finally, we examine the trilinear Higgs self-couplings. Using the results of Ref. [12] (see also 
Ref. [48]), the three-Higgs vertex Feynman rules (including the corresponding symmetry factor 
for identical particles but excluding an overall factor oii) are given by: 


9hAA ~ ”'*^[(■^3 + -^4 ~ ) 

9hAA ~ [(-^3 + -^4 ~ Z^)cg-ct ~ ZjSy-^ , 


(59) 

(60) 


9hHH — -^v[ZiSg-aC^p_^ + Z^bSy-a (5 - C^-a) + Z^Cy-ai} - 3s^y_J + Z7S^y_^Cy-a] , (61) 
9Hhh ~ ~^'^\_ZlCy-aSy_^ + Z^4^^Cy-a (3 “ '®/3-q) “ Z^Sy-^il — 3Cy_^) — Z7Cy_^Sy-a\ , (62) 


9hhh — [ZiSy_^ + Z345Sy_aCy_a + ^ZeCy^aS^^a + ^7c|-a] , (63) 

9hHH = -3^ [ZlCy-a + Z^4^Cy-^sl_^ - 3ZQSy-^cl_^ - ^T^-a] > (64) 

9hH'^H~ ~ ^ \Zj,Sy—a + Z’jCy—a^ , (65) 

9hH+H~ ~ ~'^\_Z3Cy—Q, Z-jSy—al ■ (66) 


The trilinear Higgs couplings expressed in terms of the physical Higgs masses are given in 
Appendix B. 

Consider the alignment limit, cy_a ^ 0, where h is SM-like. Then Eqs. (44) and (63) yield,® 


n — 

yhhh — yhhh 


, 2^6 f Z345 

1 + ^^Cy_a + 


2Z| 3 


Zi V 2 

where the self-coupling of the SM Higgs boson is given by 


“ o ) ^l3-a + 62(c^_«) 


9lfh = 


3ml 


(67) 


( 68 ) 


^Eq. (67) is obtained in the convention where s^_q, is non-negative, i.e. is close to 1. 
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Note that in the alignment limit, ~ [cf. Eq. (41)], which implies that Zi ~ 0.26. 
It is convenient to make use of Eq. (47) [in a convention where > 0] to write 


CjS—a ^-^6 ; 


(69) 


where 


V = 


- Ziv^) 

Inserting Eq. (69) in Eq. (67) yields 


0 ( 1 ). 


o 


m 


H 


for m\j ~ Cl(n^), 


<C 1 , in the decoupling limit. 


(70) 


Qhhh - 9hfhi 1 + 


(^345 -lZi)rf -2ri 


^2 


+ OirfZl) + Oir^^Zt) \ . 


(71) 


In the decoupling limit (where <C 1), 

ta=fi{l-^ + 0(9%“)}. 


(72) 


It follows that Qhhh is always suppressed with respect to the SM in the decoupling limit.This 
behavior is confirmed in our numerical analysis. In contrast, in the alignment limit without 
decoupling, \Zq\ is significantly smaller than 1 and ~ 0(1). It is now convenient to use 
Eq. (27) to eliminate .^ 345 , 


Qhhh — 9f^h\ i- + 


(Zy tan 2/5 — — 2 rj 


72 

Zi 


73 

(2 cot 2/5 — tan 2/5)r/^-^ + 0(Zq) 


(73) 


where the term above designated by 0(Zg) contains no potential enhancements in the limit 
of S 2 ^ — 7 - 0 or C 2/3 — 7 - 0. Given that rj ^ 0(1) in the alignment limit without decoupling, the 
form of Eq. (73) suggests two ways in which Qhhh can be enhanced with respect to the SM. Eor 
example if tan/5 ~ 1 , then one must satisfy (Zy — Zg)//tan2/5 ^ 2 + ^Zitj. Alternatively, if 
tan/5 ^ 1 , then one must satisfy Ze?/cot 2/5 ^ 1 + \Zir] (the latter inequality requires Zq < 0 , 
since cot2/5 < 0 when < /5 < In both cases, Qhhh > gfi^h possible even when \Zq\ 
and |Zy| are significantly smaller than 1. Indeed, both of the above alternatives correspond to 
Z 345 Zi and / 7 Z 345 1 in Eq. (71). 

As a second example, consider the hAA coupling given in Eq. (59) [or Eq. (B. 6 )]. Using 
Eq. (27), we find that in the alignment limit, 


9hAA = --\ml - 2 Z 5 U - (Ze - Zy)u tan2/5 + 2Z6U cot 2/5 + 0{ci3_a) 
= — TTif^ — 2 A 5 U^ + 2 Zqv'^ cot 2/5 + 0{cp—(y) [■, 


(74) 


the double decoupling limit where C 1 and \Zq\ 1, Eq. (72) shows that the deviation of Qhhh from 
the corresponding SM value is highly suppressed. 


11 











A similar computation yields the Hhh coupling given in Eq. (62) [or Eq. (B.9)], 

QHhh = - [ml - 4Z6U^cot2/5 + 2{Ze - tan2/5]c^_„ + C>(cJ_„)| . (75) 

In the alignment limit without decoupling, the 0{1) terms in Eqs. (74) and (75) that are 
proportional to Zq should be regarded as terms of 0 {ci3-a) [cf. Eqs. (69) and (70)]. That 
is, the decoupling limit [with Zq ~ 0 {1)] and the alignment limit without decoupling can be 
distinguished in the trilinear Higgs couplings. Indeed, the Hhh coupling is suppressed in the 
alignment limit without decoupling, whereas it can be of 0{v) in the decoupling limit. All 
the other trilinear Higgs self-couplings can be analyzed in the alignment limit following the 
procedure outlined above. 

Last but not least, it is noteworthy that 


9hH+H - —v[Zs + 0{c^-a)\ , (76) 

approaches a finite nonzero value in the alignment limit, with or without decoupling. This 
is relevant for the analysis of the one-loop process h —?• 77 , which has a contribution that is 
mediated by a loop. In the decoupling limit, the charged Higgs loop amplitude is suppressed 
by a factor of C>(u^/m^±) relative to the and the top quark loop contributions. But, in the 
alignment limit without decoupling, the charged Higgs loop is parametrically of the same order 
as the corresponding SM loop contributions, thereby leading to a shift of the ^ 77 decay 
rate from its SM value. This is in stark contrast to the behavior of tree-level Higgs couplings, 
which approach their SM values in the alignment limit with or without decoupling. That is, 
the loop-corrected Higgs couplings to SM particles approach their SM values in the decoupling 
limit, but can yield deviations in the alignment limit without decoupling due to internal loops 
involving light non-SM-like Higgs states. 

Before concluding this section, we examine a second theoretical distinction between the 
decoupling limit and alignment limit without decoupling. The SM Higgs sector is famously 
unnatural [50,51]. In particular, a fine tuning of the Higgs sector squared-mass parameter 
is required in order to explain the observed value of the vacuum expectation value (vev), 
V ~ 246 GeV. The 2HDM generically requires two separate and independent fine tunings. In 
addition to identifying v ~ 246 GeV, which fixes the values of li and in Eq. (31), one must 
also perform a second fine-tuning to fix the squared-mass parameter Y 2 to be of 0{v‘^). Thus, 
the regime of the decoupling limit (where Y 2 ^ v^) is less fine-tuned than the general 2HDM, 
since the natural value for Y 2 is the ultraviolet cutoff of the theory beyond which new physics 
presumably enters. As long as the heavier Higgs scalars (whose squared masses are of order 
Y 2 ) are sufficiently massive, then h will be SM-like.^^ 

In contrast, in the case of alignment without decoupling (or in the double decoupling limit), 
we have \Zq\ <C 1, which is a finely-tuned region of the 2HDM parameter space (beyond the 
two tunings discussed above) unless we can demonstrate that Ze = 0 is a consequence of an 
enhanced symmetry of the theory. The possibility of a natural implementation of alignment 

general, ^ \Ze\v'^ is sufficient to guarantee SM-like h couplings. However, in the 2HDM with 
Type-II Yukawa coupling and tand > 1, a SM-like h coupling to down-type quarks and leptons requires 
m'jj » l^elu^tand, leading to the phenomenon of delayed decoupling [12,52-54] at large tan/3. This is a 
special case of delayed alignment introduced below Eq. (58). 
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has been previously treated in [55]. In the absence of Higgs-fermion Yukawa couplings, it 
is sufficient to consider the symmetry properties of the scalar potential. Note that we have 
already imposed a softly-broken Z 2 symmetry, which yields Ae = A 7 = 0 in the original basis. 
In addition, we observe that Zq = Zt = 0 [which also implies that >3 = 0 in light of Eq. (31)] 
corresponds to an exact Z 2 symmetry in the Higgs basis. 

The conditions Zq = Z-j = 0 can be implemented in three ways. If S2/3 = 0, then only 
one of the two Higgs fields acquires a non-zero vev. This means that our original basis and 
the Higgs basis coincide (in a convention where Hi denotes the Higgs field with the non-zero 
vev), in which case the original Z 2 symmetry is unbroken. If Ae = A 7 = 0 and S 2 / 3 C 2/3 7 ^ 0, 
then setting Zq = Z-j = 0 in Eqs. (24) and (25) yields Ai = A 2 = A 345 . Such a scalar potential 
exhibits a softly-broken CP3 symmetry, one of the three possible generalized CP symmetries 
that can be imposed on the 2HDM [56].^^ Finally, if the scalar potential exhibits an exact 
CP 2 symmetry, or equivalently there is a basis in which the Z 2 discrete symmetry ('hi —)■ -|-$i, 
^2 —'^’ 2 ) and a second Z 2 interchange symmetry (•hi i—> 'h 2 ) coexist [46,56], then it follows 
that Ae = A 7 = 0, Ai = A 2 (with A 5 real), mfi = 17122 "^12 = 0- ^^is case, Eqs. (3) and 

(4) yield tan/3 = 1 .^^ The latter can be maintained when the CP 2 symmetry is softly broken 
such that m ^2 7^ 0- Using Eqs. (24) and (25) then yields Zq = Z-j = 0. Thus, in the absence 
of the Higgs-fermion Yukawa couplings, Ze = 0 is a consequence of an enhanced symmetry of 
the scalar potential, in which case the regime of alignment without decoupling and the double 
decoupling regime are both natural in the sense of’t Hooft [57]. 

If we now include the Higgs-fermion Yukawa coupling, we can still maintain the symmetry 
of the scalar potential in special cases. If the Z 2 symmetry transformation is defined in the 
Higgs basis such that H 2 is odd (i.e. H 2 —>■ —H 2 ) and Hi and all fermion and vector fields are 
even, then the resulting model corresponds a Type-I 2HDM with S 2/3 = 0, which we recognize 
as the inert 2HDM [58,59]. Indeed, if we perturb the inert 2HDM by taking Zq and Z^ small, 
then either h or H will be approximately SM-like. In the case of S 2/3 7 ^ 0, we would need to 
extend the (softly-broken) CP3 or CP2 symmetry of the scalar potential to the Higgs-fermion 
Yukawa sector. As shown in [60], no phenomenologically acceptable CP2-symmetric model 
exists. A unique softly-broken CP3-symmetric 2HDM does exist with an acceptable fermion 
mass spectrum; however this model does not appear to be phenomenologically viable due to 
insufficient CP-violation and potentially large FCNC effects [60]. Hence, for generic choices of 
the 2HDM parameters, the regime of alignment without decoupling and the double decoupling 
regime must be regarded as more finely tuned than the generic 2HDM. 

3 Setup of the numerical analysis 

In this section, we give details on the numerical procedure. In particular, we describe 
the scan of the 2HDM parameter space and the different constraints coming from theoretical 

^^If 17112 = 0 in Eq. ( 1 ) in addition to Xq = Xj = 0 , then the Z2 discrete symmetry ($1 ^ ^2 ^ —^2) 

is exact. In this case, Zq = Zj = 0 implies that Ai = A2 = A345 and mh = 17122 [^he latter via Eq. ( 30 )], and 
corresponds to an exact CPS symmetry of the scalar potential. This restriction of scalar potential parameters 
has also been obtained in [ 55 ]. 

^^Here we assume that Ai ^ A345; otherwise, the CP2 symmetry is promoted to the CPS symmetry previously 
considered. 
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requirements, signal strengths of the observed 125 GeV Higgs state, flavor physics and direct 
searches for extra Higgs states. 

Imposing a softly broken Z 2 symmetry ('hi —>■ +'hi, ^2 —$ 2 ) on the scalar potential 

given in Eq. (1) which sets Ae = A 7 = 0, the free parameters of the 2HDM scalar potential 
can be chosen to be the four physical Higgs masses the mass term mi 2 , the 

ratio of the two Higgs vacuum expectation values tan P and the mixing angle a of the CP-even 
Higgs squared-mass matrix. In this study, we choose the following ranges for the scan, 

a e [ - 7 r/ 2 , 7 r/ 2 ], tan/5 G [0.5,60], e [-(2000 GeVf, (2000 GeV)^], 

mH± e [m*, 2000 GeV], ttia G [5 GeV, 2000 GeV], (77) 

where m* is a lower bound on the charged Higgs mass originating either from the LEP direct 
searches [61] or constraints from H-physics; mainly from the Z ^ bb (Rb), e^, Auib^, B Xg'j 
and B ^ TV constraints [38-42]. In principle both h and H can have the same properties as 
the SM Higgs and thus serve as possible candidates for the observed SM-like Higgs state. In 
this paper, we consider mh = 125.5 GeV^^, taking 

tuh e [129.5 GeV ,2000 GeV], (78) 

As mentioned in Section 2 . 1 , the degenerate case rrih ~ ran is not considered in this study. 
Instead, we require a 4 GeV mass splitting between h and H in order to avoid H contamination 
of the h signal. Since we are primarily interested in the case that the electroweak gauge 
bosons acquire most of their masses from only one of the Higgs basis doublet fields, we impose 
S/ 3 -a > 0.99, which translates into |c/ 3 _q,| < 0.14. This implies that we are allowing at most a 
1% deviation from Cp = 1. This should be compared with the expected ultimate precision for 
Cv of about 2-4% at the high-luminosity LHC, and about 0 . 2 - 0 .5% at the ILC [23,28]. 

We perform a flat random scan over this parameter space using the public code 2HDMC [30] 
for a precise state-of-the-art computation of the couplings and decay widths of the various 
Higgs states. Only points satisfying stability of the scalar potential [cf. Eq. (A. 17)], coupling 
perturbativity and tree-level S-matrix unitarity are retained. We also require the S, T, and U 
Peskin-Takeuchi parameters [62] to be compatible with their corresponding values derived from 
electroweak precision observables [63]. These constraints are also checked by means of 2HDMC. 

Next we impose constraints from the non-observation of Higgs states other than the one at 
125 GeV. From the LEP direct searches for light Higgs states, we consider the cross-section 
upper limits on > Zh/H and > Ah/H from [64] and [65] respectively. For very light 

A below 9.5 GeV, the limits from Upsilon decays [ 66 ] are important, for which we follow the 
implementation in NMSSMTools 4.6.0 [67]. Moreover, we consider the limits from GMS on light 
pseudo scalars decaying into /lAii~ [ 68 ] in the mass range tua = 5.5-9 and 11.5-14 GeV, which 
are relevant in particular in Type H models. The limits from LHG searches for additional heavy 
Higgs states are also taken into account. These include the model-independent limits from the 
searches for H ZZ^ -a- 4£ from ATLAS [69] and CMS [70] and for H —>■ ZZ^*'^ —>■ 2 i2v 
from CMS [71]. However, these limits are easily evaded in our study where it is the h that 
has Cy = Sjs-a > 0.99, while HVV couplings behave as and \cis-a\ < 0.14. (This also 

^^Having performed the parameter scans before the publication of [3] which reports a central value of the 
Higgs mass of 125.09 GeV, we use 125.5 GeV as the observed Higgs mass in this analysis. 
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holds true in view of the Moriond 2015 update of the Higgs data [72].) More important are the 
limits from H,A^tt searches in gluon-fusion or associated production with a pair of b quarks 
from ATLAS [73] and CMS [74], These are particularly relevant in the large tan/5 region 
of the Type II models where a significant enhancement of the down-type fermion coupling 
to the neutral Higgs states occurs. Finally, the limits derived from the pseudoscalar search 
A —7- Zh, h ^ bb from ATLAS [75] and CMS [76] are imposed. (Limits from other searches, 
like for A Zj [77] or hh —>■ bbbb [78], have no effect on the results.) To evaluate all 
these constraints, production of the H and A via gluon-gluon fusion (ggF) and via associated 
production with a pair of bottom quarks (bbH,bbA) are computed at NNLO QCD^^ accuracy 
using SusHi-1.3.0 [32], while the vector-boson fusion (VBF) mode for the H is computed at 
NLO with VBFNLO-2.6.3 [33]. 

Signal strengths constraints coming from the precise measurements of the properties of the 
125 GeV state are taken into account by means of Lilith 1.1.2 [31]. We require each point 
of the analysis to be allowed at the 95% confidence level (CL). The CL is derived from the 
log-likelihood ratio 

A(-21nL)(7’) =-21nL(7’/2'HDM), (79) 

where L is the likelihood constructed by Lilith using up-to-date signal strength measurements, 
V represents the set of parameters of the tested point and 2HDM the best-fit point of the 
model. The Lilith database 15.04 is used for this analysis. It contains all the latest Higgs 
signal strengths measurements from ATLAS [4,80-88] and CMS [5,70,89-95] as of April 2015 
and a combined D0 and CDF result [96]. 

The updates in this version with respect to the one published in Phys. Rev. D 92 , 075004 
(2015) [arXiv:1507.00933v3] are as follows. First, we include the updated bound on the charged 
Higgs mass in Type H, mH± > 480 GeV at 95% CL [42], based on the observed rates for the 
weak radiative H-meson decay, B (instead of the previous bound of mH± > 300 GeV). 

Second, we correct a bug in the evaluation of the ATLAS and GMS limits from H,A ^ tt 
searches [73,74]: the branching ratio for A ^ tt was set to the one of a SM Higgs with the 
same mass instead of the one computed in the 2HDM. This affects the intermediate-to-large 
tan P region of the Type H models where an enhancement of the down-type fermion coupling 
occurs. In particular it eliminates part of the “opposite-sign” Cd solution in Type H. Third, 
we include the CMS constraint [97] on neutral Higgs bosons with masses between 25 GeV and 
80 GeV, produced in association with a pair of b quarks and followed by the decay into rr. 
Finally, the GMS result [79] on the search for a new heavy resonance decaying to a Z boson 
and a light resonance, followed hy Z £~^£~ and the light resonance decaying to bb or rr 
is included. Sensitive to light resonances with masses above ~ 35 GeV, this analysis puts a 
very severe constraint on gg ^ A ^ Zh with Z —?• £~^£~ and h ^ 66 in the mn = 125 GeV 
scenario [29]. Here, in the = 125 GeV scenario, it contributes to the exclusion of A with 
masses between ~ 60 GeV and ~ 150 GeV as will be shown in the next Section. 

^®The NNLO corrections for ggF are only computed for the top quark loop, as those for the bottom quark 
loop are very small. 
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Figure 1: |c/ 3 _a| versus uih in Type I (left) and Type II (right) with log^o l-^el color code. Points 
are ordered from high to low log^Q jZ^j values. 

4 Results 

4.1 Parameters 

Let us start by reviewing the relevant parameter space. Figure 1 shows the crucial relation 
between jZ^j, |c/ 3 _q,| and ttih, illustrating the different ways alignment can occur with and 
without decoupling.^® As expected, jZ^j exhibits a clear dependence on the H-h mass difference, 
see Eq. (42), and steeply drops towards zero in the limit |c^-a| =■ 0, i.e. when the h becomes 
purely SM-like. When mn is of the order of 1 TeV, one needs to be extremely close to = 1 
to have small \Zq \—for instance jZ^j ^ 10“^ requires |c/ 3 _q,| 6 x 10“® for = 1 TeV. In 

contrast, for a lighter H the departure of S/ 3 -q- from 1 can be more important—for instance the 
same |.^6| ~ 10“® value requires |c/ 3 _q,| 2 x 10 “^ for mn = 200 GeV. It is in principle always 

possible to obtain arbitrarily small values of \Zq\ if one pushes arbitrarily close to 1. For 
the purpose of the numerical analysis, we limit ourselves to |c^-a| > 10“®; we have checked that 
this captures well all features relevant for the |c/ 3 _a| —)■ 0 limit. Interestingly, as mn becomes 
larger, we observe that the decoupling limit sets a stronger upper limit on \ci 3 -a\ than the one 
set in the numerical scan (|c; 3 _q,| < 0.14). Observing a heavy uih ^ 850 GeV at the LHC would 
provide a better-than-1% indirect determination of the h-coupling to electroweak gauge bosons 
in the framework of these scenarios. 

The range of is also interesting. In principle itia can be above or below nih^H, and 
even uia < 'rnh/2 is possible and consistent with the data [27]. However, once uih is hxed, 
the allowed range of is limited (and vice versa) as illustrated in Fig. 2. We see that in 

^®In this and subsequent figures, we give 3d information on a 2d plot by means of a color code in the third 
dimension. To this end, we must chose a definite plotting order. Ordering the points from high to low values in 
the third dimension, as done for log^o \Zq\ in Fig. 1, means that the highest values are plotted first and lower 
and lower values are plotted on top of them. As a consequence, regions with low values may (partly) cover 
regions with high values. The opposite is of course true for the ordering from low to high values. To avoid 
a proliferation of plots, in each figure we show only one ordering, trying to choose the one that gives most 
information. The figures with inverted plotting order are available upon request. 
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Figure 2: mn versus in Type I (left) and Type II (right) with the color code indicating 
the value of log^g |c/ 3 -a|- Points are ordered from high to low log;^o Th® dashed lines are 

isolines of ^ 5=4 (upper line), 0 (middle line) and —4 (lower line) for \cis-a\ = 0.015 (varying 
|c; 3 _q,| from 0 to 0.14 has no visible effect on them). 


both Type I and Type II, if the scalar H is heavy and decoupled, the same is true for the 
pseudoscalar A. Conversely, if H is light, say below 600 GeV, then also A must be below 
about 800 GeV. Furthermore, it appears that for |c/ 3 _q,| < 10 “^ (or, equivalently, small \Ze\) 
iriH < itia is favored. This can be understood from Eq. (43) [or Eq. (45)]: since the 
(or sp-a) term therein is always quite small, the mass ordering between tuh and tua is 

largely determined by the sign of Z^. The value of Z 5 , in turn, is driven by A 5 [cf. Eq. (A. 13)], 
which according to our numerical analysis tends to be negative for small C/j-q. The absence of 
points over a large region of low tuh^a in Type II is in part due to the H,A^tt limits [73,74], 
which eliminate a large swath of parameter space with ~ — 1 and m-A ~ 150-350 GeV, and 
to the CMS H ZA search [79] which eliminates points down to ~ 60 GeV (with a mild 
dependence on mn)- We note that the surviving points with < 60 GeV have tan (3 <2. In 
addition, the charged Higgs mass limit mH± > 480 GeV in Type II [42] results in the elimination 
of the remaining quadrant with tuh^a ^ 400 GeV (actually up to ran ~ 430 GeV for very light 
niA). 

The interrelation between uia, ran and is illustrated in Eig. 3. The two panels show 
tuh versus with color-coding according to mH±, with the ordering going from high (blue) 
to low (red) mH± values. While the correlation of mH± with mn and tua is somewhat different 
in Type I and Type II, in both models a light charged Higgs below 500-600 GeV requires that 
the H and A also be not too heavy, with masses below about 800 GeV. When inverting the 
plotting order of mH± (not shown), we find that for any given there is a lower limit on 
tuh and for mH± ^ 1 TeV, also mH,A are of that order. In turn, when mn and are in 
the non-decoupling regime, mH± cannot be much heavier. The absence of points in the light 
mass region mH,A ^ 400 GeV in Type H (but not for Type I), already noted in the previous 
paragraph is due to the fact that in the Type H model H-physics requires ^ 480 GeV and 
at low m-A the precision electroweak T parameter constraint would be violated if ran differs 
very much from mH±- As also mentioned above, an additional band with tua ~ 150-350 GeV 
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Figure 3: mn versus itia in Type I (left) and Type II (right) with the color code indicating the 
value of mH± ■ Points are ordered from high to low m^f±. 


is cut out by the H, A ^ tt limits. We will see later that this corresponds to a large extent to 
the “opposite-sign” Cd solution with large tan/5 in Type II. 

4.2 Couplings 

The next question to address is what variations in the couplings of the 125.5 GeV state are 
still possible in the limit of approximate alignment where Cy ~ 1. In particular, recall that 
in the scan we impose > 0.99 with ruh — 125.5 GeV, without requiring however that the 
other couplings of the h be very SM-like. To answer this question, we first show in Fig. 4 the 
dependence of the reduced couplings to (up-type) fermions, see Table 2, Cp = Cp = Cp in 
Type I [Cp in Type II) on |c^_a|. The mass of the heavier scalar H is shown as a color code. 
We see that when mn is light, for only 1% deviation from unity in Cy, Cp can deviate as 
much as about 10% (20%) from unity in Type I (Type II). Inverting the plotting order of mn 
(not shown), it is interesting to note that these deviations are largest for uih ~ 700-800 GeV 
while slightly more constrained for lighter mn- On the other hand, in the decoupling limit 
the deviations in Cp are more constrained, with a maximum of 5% for mn ^ 1.2 TeV in both 
Type I and Type II. It is also interesting to observe how quickly alignment leads to SM-like 
couplings: for \cis-a\ 10“^ the deviations in Cp are limited to just a few percent no matter 
the value of mn- 

The situation is quite different for the coupling to down-type fermions, Cp, in Type II, see 
Fig. 5. First of all, the possible deviations are larger than for Cp, with Cp ranging from about 
0.70 to 1.15 even for \cis-a\ ^ 10“^. Indeed, this is an example of the delayed alignment limit 
discussed below Eq. (58); one needs |c^-a| as low as about 3 x 10“^ to have Cp within 2% of 
unity. This drives the whole phenomenology of the scenario: as we will see, sizable deviations 
of Cp from 1 lead to possible large deviations in the signal strengths even for quite small |c/ 3 _q,|. 
Inverting the plotting order of ran (not shown), we note, however, that for any given |c/ 3 _a| 
of a few times 10“^ or smaller, Cp is limited to be closer to 1 when ran is small than in the 
decoupling case with large niH- 
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Figure 4: \cis-a\ versus Cp in Type I (left) and \cis-a\ versus in Type II (right) with ran 
color code. Points are ordered from low to high mn- The points with Cp ~ 1 and \cis-a\ > 0.03 
are the points for which Cp ~ —1, i.e. the opposite-sign Yukawa coupling points, see Fig. 5. 



Figure 5: \ci 3 -a\ versus Cp in Type II with uih color code for the full Cp range (left) and 
zooming on the Cp > 0 region (right). Points are ordered from low to high mn- 


Moreover, Cp = 1 is not possible unless |c; 3 _a| is very small (again a few times 10“^ or 
smaller) as a consequence of the lower bound tjs > 0.5 imposed in the analysis. Large pos¬ 
itive deviations of Cp, up to ~ 1.12, would indicate rriH ^ 750 GeV. On the contrary, Cp 
values which are substantially smaller than 1 can be achieved in both the decoupling and non¬ 
decoupling regimes except for a small island of points located around Cp ~ 0.8 and |c/ 3 _q,| 0.1 

that is achieved only for ran ^ 400 GeV. Thus, for instance, a discovery of a light H state in 
association with a measured value of Cp ~ 0.8 would give an indirect way to probe sub-percent 
deviation of Cy in this Type II scenario. 

Finally, for light tuh the sign of Cp relative to Cy and Cp can be opposite to the cor¬ 
responding SM value. This is realized for not so small values of |c/ 3 _a| > 0.07, i.e. at the 
boundary of what we consider as the alignment limit, for 330 GeV < ran < 660 GeV, 
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Figure 6: Fermionic couplings versus tan P in Type I (upper panel) and Type II (lower panels) 
with \ci 3 -a\ color code. Points are ordered from high to low |c^_a|. 


350 GeV < rriA < 660 GeV and 0.22 <\Zq\< 0.90. For the points in this region, the up-type 
coupling is very close to 1, corresponding to the few isolated points observed in the right panel 
of Fig. 4. As discussed in [54], the eventual LHG Run 2 precision will allow one to either con¬ 
firm or eliminate the opposite-sign coupling possibility using precise signal rate measurements 
of the h in a few channels. Should the opposite-sign coupling be confirmed, one would expect 
to also see A signals (plus perhaps H signals) in the above mass range, thereby providing a 
confirmation of this scenario. (The cross sections for A and H signals will be discussed in 
Section 4.4.) It should also be noted here that this region is much reduced by the correction of 
the H,A^ TT constraints as compared to the previous version [arXiv:1507.00933v3]. 

The tan /3 dependence of the fermion couplings of h is shown in Fig. 6. We see that large 
tan P leads to very close to 1 in Type I and very close to 1 in Type II. However in Type II, 
at large tan/5, small is not enough to drive —?• 1: the approach to SM-like coupling 

is delayed, as discussed in Section 2 in the text below Table 2. Note also that the opposite- 
sign solution in Type II requires tan/5 > 10 and Cp ~ 0.9994 (which is experimentally 
indistinguishable from exact alignment). 


20 


















Figure 7: |c^-a| versus in Type I (left) and Type II (right) with tuh color code. Points are 
ordered from low to high mn- 


The loop-induced coupling to photons, is presented in Fig. 7. Even at very small c^_a, 
C!^ can deviate substantially from 1. This is due to the charged-Higgs contribution to the hyy 
coupling. This contribution can be large with either sign, positive or negative, in Type I, while 
in Type II large contributions are always negative and suppress [54] . Note in particular the 
Type II points with ~ 0.95 associated with the opposite-sign (Fy, cases for which the charged 
Higgs loop contribution does not decouple and always leads to a suppression. Regarding the 
loop-induced coupling to gluons, in the Type I model, is equal to Cp (up to NLO), the 
dependence of which on {cg-al was presented in Fig. 4. In the case of Type II, Cj and Cl^ 
are very similar despite the difference between up and down-type couplings, this being due 
to the fact that the 6-loop contribution to is rather small. The one exception in the case 
of Type II arises for the opposite-sign scenario for which the 6-loop contribution changes sign 
and interferes constructively with the t-loop contribution. In this case, Cg is always enhanced, 
Cj - 1.06 [54]. 

While the exceedingly small deviations in Cy that we consider here will most likely not be 
directly accessible at the LHC, precision measurements of the other couplings together with 
a measurement of, or a limit on, itih^a can be used for consistency checks and for eventually 
pinning down the model. Of special interest in this context is also the triple Higgs coupling. 
The dependence of Chhh = 9hhh/gf^h and mn is shown in Fig. 8. It is quite striking 

that large values of Chhh > 1 (up to Chhh ~ 1-7 in Type I and up to Chhh ~ L35 in Type H) 
can be achieved in the non-decoupling regime, roughly mn ^ 600 GeV, for Icg-d values of 
the order of 0.1, whereas for heavier uih, Chhh is always suppressed as compared to its SM 
prediction. These features were explained in the discussion below Eq. (67).^^ Note also that 
for rriH ~ 1 TeV, Chhh approaches the SM limit of 1 as \cis-a\ decreases more slowly than is the 
case for lighter substantial deviations Chhh < 1 are possible as long as |c/ 3 _q,| is roughly 
greater than a few times 10“^. This comes from the (2ZqIZ i)c^-a, term in Eq. (67): since, in 
the convention where sg-a > 0, Z^cg-a is always negative, cf. Eq. (46), and since Zq can be 


^^This cannot be seen directly in Fig. 8, but we verified that points with niH > 630 GeV never have Chhh > 1- 
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Figure 8: \cis-a\ versus the reduced triple Higgs coupling Chhh in Type I (left) and Type II 
(right) with tuh color code. Points are ordered from high to low mn values. 


2HDM Type I, =125.5 GeV 2HDM Type II, =125.5 GeV 



9HHh [GeV] [GeV] 


Figure 9: |c; 3 _a| versus the triple Higgs coupling gnhh in Type I (left) and Type H (right) with 
niH color code. Points are ordered from high to low mn values. 


sizable when uih ~ 1 TeV, see Fig. 1, this can lead to a suppression as extreme as Chhh ~ 0.1. 
(For mn ^ 1 TeV the deviations are smaller in part because the possible range of Cjs-a is 
limited as seen in Fig. 1.) For very light mn, on the other hand, Zq is much smaller and hence 
the deviations with Chhh < 1 are more limited. For ran ^ 250 GeV we find Chhh ~ 0.80-1.40 
in Type I and Chhh ~ 0.95-1.13 in Type H. This is at the limit of what can be measured, as 
the expected precision is about 50% at the high-luminosity options of the LHC and the ILC 
with 500 GeV, and about 10-20% at a 1-3 TeV e+e“ linear collider with polarized beams [28] . 

The relation between the triple Higgs coupling gHhhi \ci3-a\ and mn is presented in Fig. 9. 
In Type I, large values of gHhh can be achieved in the non-decoupling regime for \cis-a\ of the 
order 10“h This is also true in Type H, though the range of gHhh is somewhat smaller. We 
observe moreover that for given |c/ 3 _a| C 10“^, the achievable Hhh coupling grows with uih- 
Nonetheless, as will be shown in Section 4.4, the H ^ hh decay is mostly relevant below the 
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tt threshold. Moreover, in the exact alignment limit, the Hhh coupling vanishes. 

4.3 Signal strengths 

The variations in the couplings to fermions discussed above have direct consequences for 
the signal strengths of the SM-like Higgs boson. Since the results depend a lot on the fermion 
coupling structure, we examine this separately for Type I and Type II. 

Let us start with Type I. Figure 10 shows the signal strengths for gluon-gluon fusion and 
decay into 77 left panel), and decay into ZZ* {iigg{ZZ*), right panel). Recalling 

that Cp varies between 0.87 and 1.11 in Type I and comparing with Fig. 7, it is clear that 
the variation in iJ^ggi'yj) comes to a large extent from the charged Higgs contribution to the 77 
loop. Even for |c/ 3 _a| —?• 0, large deviations from 1 can occur due to a sizable charged Higgs 
contribution or the presence of a light pseudoscalar uia < mh/‘i that increases the SM-like 
Higgs total width. On the other hand, in the decoupling limit, the charged Higgs loop is small 
and C!^ is largely determined by the relative size of the top and bottom loops compared to the 
W loop (which enters with opposite sign). On the contrary, Cg is solely determined by the 
size of the t and b loop contributions. One finds numerically that the hyy coupling is more 
suppressed than the hgg coupling is enhanced, so that ^ 1 in the decoupling regime. 

In contrast, ij!^g{ZZ*) shows less variation, gJ^g^ZZ*) = [0.92,1.04] if the h —>■ AA decay 
channel is closed, with small excursions around 1 allowed in the decoupling limit. It also exhibits 
a less distinct dependence on mn compared to gggi'j'j). The reason is that g!^g{ZZ*) is driven 
by Cp and tan/5, as illustrated in Fig. 11 . The dependence on Cp is clear as larger (smaller) 
Cp leads to larger (smaller) cross section for gg —>■ h. The dependence on tan (3 results from an 
interplay between the top (which drives the gg ^ h cross section) and bottom (which drives 
the total h width) Yukawa couplings both given by Cp = The scattered points 

with suppressed g!^g{ZZ*) are those where the h AA decay mode is open and increases the 
total width. An analogous picture emerges for the VBF-induced hrr signal strengths, since 
= b^^ggiZZ*) in Type I. 

In Type H, we find that the situation is quite different. Here, the signal strengths are driven 
by both the top quark coupling, which impacts Cg, and by the bottom Yukawa coupling Cp, 
which also enters Cg and, often of greatest importance, determines the h ^ bb decay width. 
In Fig. 12 we therefore show the signal strengths gggi'j'j), g!^g{ZZ*) and Type H 

comparing the dependence on mn (left panels) to the dependence on \Cp\ (right panels). Note 
that the ran dependence of the signal strengths reflects the ran dependence of Cp in Fig. 5. As 
a consequence, and g^g{ZZ*) can be enhanced in the decoupling regime, with values 

going as high as 1.4-1.5 (mainly due to suppression of the total h width), to be compared to the 
current model-independent 95% CL limits of gggi'j'j) G [0.76,1.69] and g^g{ZZ*) G [0.71,1.80]. 
Suppression is also possible, reaching a level of IJ^ggi'J'y) =0.74-0.76 for low mn if |c/ 3 -a| > 0.01 
but limited to 0.9 for large uih ^ 1250 GeV. For all mn, the amount of possible suppression 
decreases systematically with decreasing |c/ 3 _q,|. For /Uvbf('^'^) behaviour is exactly opposite. 
For completeness we note that the horizontal bar at |c; 3 _a| ~ 10“^ is the Cp < 0 region, and 
the scattered points are those where the h —?• AA decay is open.^® Finally note that as |c^-a| 
decreases, the signal strengths in Type H converge to 1 much more slowly than in Type 1. 

^®This region has been sizably affected by the updated constraints as described at the end of Section 3. 
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Figure 10: Signal strengths in Type I for the 125.5 GeV state, for gg ^ h ^ (left) and 
gg ^ h ^ ZZ* (right) with mn color code. Points are ordered from low to high mn values. 
Points with figg(ZZ*) <0.92 are ones for which h —?• AA decays are present, so that the total 
h width is increased, which suppresses this particular channel’s rate. 



Figure 11: Signal strength for gg ^ h ^ ZZ* in Type I for the 125.5 GeV state with Cp (left) 
and tan P (right) color code. Points are ordered from low to high Cp and tan P values. 


This is a consequence of the delayed alignment of to 1 in Type II when tan P is large. An 
additional effect arises in due to the charged Higgs loop contribution to the h —?• 77 

amplitude. In particular, there exists an intermediate range of charged Higgs masses^® for which 
QhH+H- ~ —2m^±/u [cf. Eq. (B.12)], which yields a constant non-decoupling contribution that 
suppresses the /i —>■ 77 amplitude [54] (see also [98,99]). Indeed, even for values of \ci 3 -a\ as low 
as 10 “^, this signal strength does not converge to 1 until uih (and thus mH±) is above about 
1 TeV. 

Putting everything together we find quite distinct correlations of signal strengths in both 

this intermediate mass region, the charged Higgs mass is given by Eq. (33), where Y 2 ~ 0{v^) and .^ 3^1 
such that the upper bound of Z 3 is constrained by its unitarity bound. 
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Figure 12: Signal strengths in Type II for the 125.5 GeV state with mn (left) and \C^\ color 
code. Points are ordered from low to high mn and \C^\ values. 
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Figure 13: Correlations of signal strengths in Type I, on the left illustrating the dependence 
on niH, on the right illustrating the dependence in |c/ 3 _a|. Points are ordered from low to high 
niH values (left) and high to low |c/ 3 _q,| values (right). 


Type I and Type II that depend on whether the additional Higgs states are decoupled or 
not. This is illustrated in Fig. 13 for Type I and in Fig. 14 for Type II. In both figures, 
the panels on the left show the dependence on ttih while the panels on the right show the 
dependence on |c/ 3 _q,| for the non-decoupling regime with mn < 600 GeV. We note that there 
are definite combinations of signal strengths that cannot be reached in the decoupling regime. 
A measurement of such values would be a very strong motivation to look for additional light 
Higgs states. In turn, when the masses of additional light Higgs states are measured, signal 
strength correlations as shown in Figs. 13 and 14 can help pin down the model. Furthermore, 
for niH < 600 GeV even in the apparent alignment limit \ci3-a\ 0 there can be deviations in 

the signal strengths from unity that cannot be mimicked by decoupling. 

Examples for Type I are the suppression of both hgg(77) ^de combination 

fJ-ggill) > 1 with figg{ZZ*) 1. The former case is also present in Type H for light uih, 
while the latter does not occur at all in Type H. More concretely, in the decoupling regime 
of Type H, ~ whereas for fight uih one can have < ii^g{ZZ*) 
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Figure 14: Correlations of signal strengths in Type II, on the left illustrating the dependence 
on niH, on the right illustrating the dependence in |c/ 3 _a|. Points are ordered from low to high 
niH values (left) and high to low |c/ 3 _q,| values (right). 
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even if |c^-a| is very small (comparing Fig 14, top row, left vs. right). Another example 
is the simultaneous suppression or enhancement of and Pvbf( 77) Type I, that is 

not possible in the decoupling regime (cf. Fig 13, bottom left). In Type II, one can have 
a simultaneous enhancement, up to 1.45 of and Pvbf( 77) decoupling regime, 

but simultaneous suppression is limited to ~ 0.9-0.95 (cf. Fig 14, middle left); simultaneous 
suppression to a level of ~ 0.8 is however possible in the alignment limit for mn ^ 300 GeV, 
i.e. well away from the decoupling regime. Precise enough signal strength measurements could 
therefore provide strong hints that we are in the alignment without decoupling regime of a 
2HDM even if no additional Higgs states are discovered at that time. 

4.4 Cross sections for H and A production 

Let us now turn to the prospects of discovering the additional neutral states. The two 
largest production modes at the LHC are gluon fusion, gg ^ X, and the associated production 
with a pair of 6-quarks, bbX, with X = A, H. The correlations of the gg ^ X and bbX cross 
sections at the 13 TeV LHC in the non-decoupling regime mn < 600 GeV are shown in Fig. 15 
for the Type I model and in Fig. 16 for the Type H model. We show the points that pass all 
present constraints (in beige) and highlight those that have a very SM-like 125 GeV Higgs state 
by constraining all the following signal strengths to be within 5% or 2% of their SM values, 
respectively, denoted as SM±5% (red) and SM±2% (dark red): 

)) f^ggi'^A^ lA^ BFill) ^ )) iAbfAAi hyir(^^)) ' (®0) 

We start the discussion with production of A in Type I, shown in the left panel of Fig. 15. 
There is a strong correlation between the two production modes, gluon fusion and 66 associated 
production, which stems from the fact that the relevant couplings are the same up to a sign: 

— —C^ = cot /3. The larger spread in a{bbA) observed for a{gg —>■ A) > 10“^ pb comes 
from the fact that for < 2mt GeV the bbA cross section grows faster with decreasing 
than that of gg —>■ A. Therefore, along a line of fixed <7{gg A) in the plot, a point with 
higher a (bbA) has a smaller m^. Note also that there is an interference of the top and bottom 
loop diagrams in gg ^ A which changes sign depending on m^. Overall, however, a{gg —?• A) 
is always at least two orders of magnitude larger than a {bbA). 

The points with largest cross sections, a{bbA) 10 pb and a{gg ^ A) a; 1000 pb, corre¬ 
spond to the case uia < rnh/^ which was studied in detail in [27]. One feature of this region is 
that lAgill) lAg{ZZ*,WW*) always differ from each other by about 10%. Gonstraining 
all h signal strengths of Eq. (80) within 5% of unity therefore eliminates these points. Other 
points with high cross sections, but not in the very light pseudoscalar region, would also be 
eliminated by the SM±5% or SM±2% requirements. However, in this non-decoupling regime 
of niH < 600 GeV, points with sizeable cross sections up to 0.2 pb for a{bbA) and up to about 
40 pb for a{gg —?• A) still remain even at the SM±2% level. At this same SM±2% level, the 
smallest a{gg —?• A) is about 0.1 fb. 

Regarding production of the scalar H in Type I, shown in the right panel of Fig. 15, the 
correlation is even stronger between a{bbH) and <j{gg —)■ H) since both are driven by the same 
fermionic coupling Cp = sin a/ sin /5. Note that, as in the A case, the gluon-fusion cross section 
is always larger than that for 66 associated production. Sizable cross sections are still allowed 
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Figure 15: a{hbX) versus a{gg —?• X) for X = ^ (left) and X = H (right) in Type I at the 
13 TeV LHC for points satisfying all present constraints (in beige) as well as points for which 
the signals strengths from Eq. (80) are within 5% and 2% of the SM predictions (in red and 
dark red, respectively). The dashed lines indicate a{bbX) = a{gg —>■ X). 
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Figure 16: As Fig. 15 but for Type II. 


under the SM±2% constraint, which implies that in the non-decoupling regime there is a strong 
possibility of detecting a non-SM-like scalar state at the LHC. The structure of Cp is however 
such that the coupling can equally well be very much suppressed, leading to extremely small 
cross sections. We will come back to this below. 

The corresponding results for Type H are presented in Fig. 16. In contrast to Type I, both 
bb associated production and gluon-gluon fusion modes for Type H are in principle important 
since either can be dominant in different regions of the parameter space. There is only modest 
correlation between the two production modes due to the more complex structure of the Type H 
fermionic couplings. For A production, one clearly sees the < mh/2 region as the detached 
scattered points with very large cross sections. As for Type I, these points disappear under the 
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Figure 17: Cross sections in Type I (left) and Type II (right) for gg ^ X at the 13 TeV LHC 
as functions of mx for X = H (upper row) and X = H (lower row) with tan /3 color code. In 
all four plots, points are ordered from low to high tan /5. 


SM±5% constraint. Still, even for SM±2%, cross sections for bbA close to 1 pb and for gg ^ A 
of a few pb can be achieved (although not simultaneously). For H production a similar picture 
emerges, with the maximal cross sections however being a factor of a few smaller than for A 
production. The minimal cross sections in this tuh < 600 GeV non-decoupling regime for the A 
and H are correlated in a way that is very favorable for discovery during Run 2 of the LHC. For 
example, if cr^gg A) takes on its minimum SM±2% value of a few fb then a{bbA) > 20 fb, 
whereas if a{bbA) takes on its minimal value of fewxlO“^ fb then a{gg ^ H) 300 fb. These 
cross section levels imply that the A should be discoverable in at least one of the two production 
modes even in the extreme alignment limit. Analogous arguments hold for H production. 

Before considering specific decay channels of A and H, we present in Fig. 17 the gluon- 
fusion cross sections in Type I and Type H as functions of ttia and mn at the 13 TeV LHC. 
Here, the color code shows the dependence on tan In Type I, the gg ^ A cross section is 

avoid a proliferation of plots, we choose to show here only the results for gluon fusion; all corresponding 
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proportional to cot^/3; this explains why it is larger (smaller) at lower (higher) tan/3. A cross 
section of 1 ( 0 . 1 ) fb is guaranteed for tua as large as ~ 600 (850) GeV. On the other hand, 
the gg ^ H cross section in Type I is proportional to {C^Y and can take on extremely small 
values for uih ^ 850 GeV. The reason is that, in this region, the reachable values of C/j-q, are 
high enough such that a cancellation between the two terms of = (s^-a — c^-a/^/?) occurs 
and leads to an almost vanishing coupling. In contrast, for tuh ^ 850 GeV, this cancellation 
is not possible as the values of Cjs-a are forced to be smaller as can be seen in Fig. 1. In 
Type II, the A production cross section can be very large in the very low tua region as noted 
in [27] and any mass smaller than 1.1 TeV gives a. gg ^ A cross section larger than 1 fb. For 
gg ^ H, a cross section > 1 (0.1) fb is guaranteed up to mn ~ 850 GeV ( 1.2 TeV). From 
these considerations the prospects for discovering the additional neutral states look promising 
should alignment without decoupling be realized. 

Let us now turn to specihc signatures. Figure 18 presents the cross sections ior gg ^ A ^ Y 
for Y = 77 , rr, tt in Types I and II. Note that the y-axis is cut off at 10“^ pb. Although much 
lower values of the cross section are possible, we do not show these lower values since they will 
certainly not be observable at the LHC. As expected, for the 77 and rr final states, the cross 
sections fall sharply above the ti threshold, with the noticeable exception of the A ^ tt decay 
in Type II due to the strong constraints from LHC direct searches that exclude points with 
large corresponding cross section. For the A ^ 77 decay, cross sections of 0.1 fb are reachable 
for niA A 470 GeV (m^ < 530 GeV) in Type I (H) but not guaranteed. The maximal cross 
section is ~ 30 fb in Type I and ~ 100 fb in Type H (not considering the < mh/2 region). 
In both Types I and H, the gg ^ A ^ tt cross section can be substantially larger. In Type I, 
0.1 fb is reachable for ttia Y, 600 GeV, while in Type H tua Y: 550 GeV guarantees a cross 
section larger than 0.1 fb. In both cases, very large cross sections are predicted at low ttia- The 
gg ^ A ^ tt cross section peaks around 100 pb in both Types I and H and is guaranteed to be 
larger than 0.1 fb in Type H for 350 < < 600 GeV. These sizable cross sections therefore 

provide interesting probes of the extended Higgs sector in the alignment limit. 

The corresponding results for the H cross sections are presented in Fig. 19. Sizable values 
of a X BR are possible in both Types I and H but heavily suppressed values are still possible for 
most of the cases. Only in Type H, for H ^ tt (as well as for —)■ ft), is the corresponding 
cross section guaranteed to be larger than 0.1 fb for mn ^ 460 GeV {itia ~ 400 GeV). Note that, 
for both Types I and H, the cross sections for A/H decays into a muon pair are related to the 
A/H — 7 - TT ones through B{A/H /r/i) ~ (m^/ruT-)^ x B{A/H tt) ~ B{A/H tt)/280. 

Non-standard production modes of the SM-like state, through A —>■ Zh and H —>■ hh, are 
presented in Fig. 20. While these can be interesting discovery modes for the A and/or H, their 
cross sections can also be extremely suppressed. For gg ^ A ^ Zh^ the tan/3 dependence, 
which follows the dependence of the gg ^ A cross section shown in Fig. 17, explains a part 
of this suppression. Moreover, the AZh coupling is proportional to which is suppressed 
in the alignment region. Nevertheless, the gg ^ A ^ Zh cross section can be > 100 fb for 
niA ^ 600 GeV in both Types I and H. The gg ^ H ^ hh cross section, as expected, attains 
its maximum below the tt threshold in both Types I and H and can reach about 10 pb at low 
tan/3. For any mn, the cross section can however also be extremely suppressed. 

A comment is in order here on the possible “feed down” (FD) [13, 100] to the production 

results for the bb cross section can be provided upon request. 
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Figure 18: Cross sections times branching ratio in Type I (left) and in Type II (right) for 
gg ^ A ^ Y at the 13 TeV LHC as functions of for T = 77 (npper panels), Y = tt 
(middle panels) and Y = tt (lower panels) with tan /3 color code. Points are ordered from low 
to high tan /?. 
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Figure 19: Cross section times branching ratio in Type I (left) and in Type II (right) for 
gg ^ H ^ Y at the 13 TeV LHC as fnnctions of mn for T = 77 (upper panels), Y = tt 
(middle panels) and Y = tt (lower panels) with tan /3 color code. Points are ordered from low 
to high tan /?. 
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Figure 20: Cross sections times branching ratio in Type I (left) and in Type II (right) for 
gg ^ X ^ Y ai the 13 TeV LHC as functions of mx for X,Y = A,Zh (upper panel) and 
X,Y = H, hh (lower panel) with tan /3 color code. Points are ordered from low to high tan {3. 


of the 125 GeV state through the decay of heavier Higgs states, which might distort the Higgs 
signal strengths. This issue was approximately addressed in section HI.C of [13] by imposing the 
“FDOK” conditions /i|^ < 0.3 and jJ^Yh+bbh < which limit the FD contamination of Zh 
associated production and of ggF+66h production to 30% and to 10% respectively, at the cross 
section times branching ratio level. Imposing these conditions here would remove the points 
with ai^{gg —>■ %) x BR(74 —>■ Zh) > 0.2 pb and (7is{gg H) x BR(iL ^ hh) > 2 pb in Fig. 20. 
This is, however, a maximally conservative constraint for two reasons. Firstly, the amount of 
FD is computed without accounting for any reduced acceptance of such events into the 125 GeV 
signal as a result of the experimental cuts used to define the gg ^ h, bbh ov Z* ^ Zh channels. 
Secondly, it puts individual limits on specific production x decay modes instead of including all 
signal strengths in a global fit, which is the approach followed in this paper. Indeed, when 
directly adding the contribution of gg ^ A ^ Zh to the Zh signal strength in the global fit, 
it turns out that only cross sections of ai 3 {gg ^ %) x BR(yl —)■ Zh) > 1.6 pb are definitely 
excluded. This still assumes that the signal acceptance of the experimental analysis is the 
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same for gg ^ A ^ Zh as for gg ^ Z* ^ Zh, which should however be a reasonable 
approximation, as the main difference would be the Zh invariant-mass distribution, which is 
not used as a selection criterion in this case. The contribution oi H ^ hh to the h signal 
strengths is a more difficult question, as here the acceptances (in each final state considered 
in the experimental analyses) will certainly be different from those of single h production. A 
detailed study based on event simulation would be necessary to better understand the impact 
of FD on the 125 GeV Higgs signal, but this is beyond the scope of this paper. 

Finally, if the mass splitting is large enough, A —)■ ZH, H —>■ ZA, and H AA decays offer 
intriguing possibilities for discovering the extra non-SM-like neutral Higgs states in the regime 
of approximate alignment without decoupling. In Fig. 21, the cross sections for gg ^ A ^ ZH, 
gg ^ H ^ ZA and gg ^ H ^ AA are exhibited. Large gg ^ A ^ ZH cross sections are 
obtained for large mA — ran splitting.Looking back at Fig. 2 one sees that, in both Type I and 
Type H, the splitting can only be large for tua ^ 650 GeV. This explains the preponderance 
of low ran points with cross sections up to 20 pb (6 pb) in Type I (H) for tua ^ 650 GeV. (In 
Type H the mH± > 480 GeV constraint allows a large enough tua — rnn mass splitting only 
for tua ^ 350 GeV.) However, gg ^ A ^ ZH can also be heavily suppressed; since the AHZ 
coupling is proportional to S/j-a, this suppression is a purely kinematical effect. 

Turning to the H —>■ ZA and H AA signatures, in Type I we observe a depleted area 
for niH > 300 GeV and cross sections of the order of 0.1 pb. In this region, tan/3 = 2-10 and 
Z^ is small or negative leading to H and A masses for which the H —>■ ZA, AA decays are 
kinematically forbidden [cf. Eq. (45)]. In the region below, tan/3 > 10 and Z 5 can be large 
enough to achieve tuh > tua + mz and/or uih > 2mA, but nevertheless the cross section is 
small because of the tan/3 dependence of a{gg —>■ H), see Fig. 17. The distinct branch with 
gg ^ H ^ ZA and gg ^ H ^ A A cross sections larger than about 1 pb, on the other hand, 
has tan/3 < 2 and A 5 ~ 0. Here, the term proportional to sin2/3 in Eq. (23) gives a large 
enough Z 5 > 0 so that the H ^ ZA and/or H —?• AA decay is kinematically allowed. The 
small tan/3 leads to a large gg ^ H production cross section, see again Eig. 17. In Type H, 
gg ^ H ^ ZA and gg ^ H ^ AA cross sections can also be large (even above 1 pb for 
H ^ ZA) in the non-decoupling regime. However, due to the charged Higgs mass constraint 
these processes are allowed only for mn ^ 430 GeV. A detailed phenomenological analysis of 
the A ZH and H ZA decays at the LHC was performed in [ 102 ]. 

Last but not least, note that due to the kinematic constraint mn > 2mA and the non trivial 
correlation between mn and observed in Eig. 2 , the H AA channel is only open for 
niH ^ 700 GeV. In Type I the branch of points with cross sections ranging from about 10“^ pb 
to 10 pb is mainly populated by < 100 GeV points with relatively low tan/3 < 10. In 
Type H, due to the A —?• rr and H —?• ZA experimental constraints, only points with low 
mA ^ 60 GeV and tan/3 < 2 experience H AA decays. All in all, this channel offers a 
complementary probe to the low region discussed in [27]. 


large splitting — niH ~ v can be motivated by the possibility of a strong first order phase transition 
in 2HDMs [101]. 
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Figure 21: Cross sections times branching ratio in Type I (left) and in Type II (right) for Higgs- 
to-Higgs signatures at the 13 TeV LHC, in the npper panel gg ^ A ^ ZH with mn color 
code, and in the middle and lower panels for gg ^ H ^ ZA and gg ^ H ^ AA, respectively, 
with m-A color code. Points are ordered from high to low itia or mn, with the exception of the 
H —7- ZA, and H —?• AA plots in Type II, where points are ordered from low to high 
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5 Conclusions 


While the Higgs measurements at Run 1 show no deviations from the SM, conceptually 
there is no reason why the Higgs sector should be minimal. Indeed a non-minimal Higgs sector 
is theoretically very attractive and, if confirmed, would shine a new light on the mechanism of 
electroweak symmetry breaking dynamics. 

In this paper we focused on CP-conserving 2HDMs of Type I and Type H, investigating the 
special situation that arises when one of the Higgs mass eigenstates is approximately aligned 
with the direction of the scalar field vacuum expectation values. In this case, the IT^ and Z 
gauge bosons dominantly acquire their masses from only one Higgs doublet of the Higgs basis. 
Moreover, the coupling of that CP-even Higgs boson to the gauge bosons tends towards the 
SM value, Cy —t 1- While this is automatically the case in the decoupling limit when the 
extra non-SM Higgs states are very heavy, such an alignment can also occur when the extra 
Higgs states are light, below about 600 GeV. We specifically investigated the phenomenological 
consequences of alignment without decoupling and contrasted them to the decoupling case. Two 
aspects are interesting in this respect: one being precision measurements of the couplings and 
signal strengths of the SM-like Higgs boson at 125 GeV, the other being the ways to discover 
the additional Higgs states of the 2HDM when they are light. 

In addition to an in-depth theoretical discussion, we performed a detailed numerical analysis 
for the case that the SM-like state observed at 125 GeV is the lighter of the two CP-even Higgs 
bosons of the 2HDM, h. In this study we allowed for 1% deviation from unity in Cy, which 
corresponds to the ultimate expected LHC precision at high luminosity. The results can be 
summarized as follows: 

1. In the alignment limit without decoupling, despite Cy being very close to 1, the fermionic 
couplings of the 125 GeV Higgs can deviate substantially from the SM values. Concretely, 

can deviate as much as about 10% (20%) from unity in Type I (Type H), and as 
much as 30% in Type H. 

2. While rather quickly approaches 1 with increasing uih and/or C/s-a 0, the approach 

of the bottom Yukawa coupling to its SM value in the alignment limit is delayed in Type H, 
with Gp 0.70-1.15 even for |c/ 3 _a| 10 Large values of > 1 are associated with 

light H, A. Moreover, for 330 GeV < uih < 660 GeV and 350 GeV < uia < 660 GeV, 
there is an allowed region with G/, —1±0.2; this “opposite-sign” solution can be tested 

decisively at Run 2. 

3. The trilinear hhh coupling can also exhibit large deviations. Large values of Cuhh > 1 
(up to Chhh ~ 1-7 in Type I and up to Chhh ~ 1-35 in Type H) can be achieved in the 
non-decoupling regime mn ^ 600 GeV, for |c^-a| of the order of 0.1, whereas for heavier 

Chhh is always suppressed as compared to its SM prediction. The suppression can 
be about 50% for uih of ~ 1 TeV and much larger for lighter rriH- 

4. For the ratios jJxiY) of the X ^ h ^ Y signal rates relative to the SM prediction, 
we found distinct correlations of these signal strengths in both Type I and Type H that 
depend on whether the additional Higgs states are decoupled or not. In fact, in the regime 
of alignment without decoupling, there are characteristic combinations of the tjfxiY) signal 
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strengths that cannot be mimicked by the decoupling limit. However, it is of course also 
possible that all signal strengths converge to 1 even though the additional Higgs states 
are very light. 

5. A decisive test of the alignment without decoupling scenario would of course be the 
observation of the additional Higgs states of the 2HDM in the mass range below about 
600 GeV. We delineated the many possibilities for such observations. While there are 
no guarantees in the case of the Type I model, in the Type H model there is always a 
definite lower bound on the gg ^ A, H ^ tt cross sections at the LHC at any given 
niA- For low tan/5 ~ 1, this lower bound is still of order 0.1 fb for ~ 500 GeV, a 
level that we deem likely to be observable at the LHC during Run 2. For high tan /5, the 
lower bound is roughly two orders of magnitude higher and only falls below the 0.1 fb 
level for mA,H ^ 1-2 TeV, which is already in the decoupling region. Moreover, while in 
Type I gluon-gluon fusion is always dominant for H ox A production, in Type H both 
bb associated production and gluon-gluon fusion modes are in principle important since 
either can be dominant in different regions of the parameter space. 

6 . Higgs-to-Higgs decays of the non-SM-like states (A —?• ZH, H ZA, H AA) also 
open intriguing possibilities for testing the regime of alignment without decoupling, with 
cross sections often in the range of 1-10 pb (although they can also be quite suppressed). 
Particularly promising are gg ^ H ^ ZA and gg ^ H ^ AA in Type H for light 
pseudoscalars below about 100 GeV; for such a light A, mn can be at most ~ 650 GeV, 
and a X B values for these channels typically range from 10 fb to 10 pb. 

In short, it is possible that the observed 125 GeV Higgs boson appears SM-like due to 
the alignment limit of a multi-doublet Higgs sector. The alignment limit does not necessarily 
imply that the additional Higgs states of the model are heavy. Indeed, they can be light and 
non-decoupled and thus lead to exciting new effects to be probed at Run 2 of the LHC. 


Acknowledgements 

J.B. and S.K. thank Jose Miguel No and Ken Mimasu for discussions on the A —)■ ZH decay 
during the “Physics at TeV Colliders” workshop in Les Houches 10-19 June 2015. 

This work was supported in part by the Research Executive Agency (REA) of the Euro¬ 
pean Union under the Grant Agreement PITN-GA2012-316704 (HiggsTools) and by the “In- 
vestissements d’avenir, Labex ENIGMASS”. H.E.H. is supported in part by U.S. Department 
of Energy grant DE-EG02-04ER41286. J.E.G. and Y.J. are supported in part by the US DOE 
grant DE-SC-000999. Y.J. also acknowledges support by the LHC-TI fellowship US NSF grant 
PHY-0969510, and he thanks the LPSC Grenoble for its hospitality. J.E.G., H.E.H. and S.K. 
are grateful for the hospitality and the inspiring working atmosphere of the Aspen Center for 
Physics, supported by the National Science Foundation Grant No. PHY-1066293, where this 
project was initiated. 


38 



A Scalar potential quartic coefficients in the Z 2 -basis in 
terms of Higgs basis coefficients 

In Eqs. (21)-(25), we have provided expressions for the Higgs basis quantities Zj in terms 
of the quartic coefficients of the scalar potential A, defined in Eq. (1). In this appendix, we 
provide the inverse of Eqs. (21)-(25) by expressing the Aj in terms of the Zj. 

Ai = Zic^ + Z2S^ + IZ345S2/3 — 2 s 2 / 3 (c ^-^6 + ; (A-1) 

A2 = ZiS^ + Z2C^ + IZ345S2/3 + 2 S 2 / 3 (-S ^-^6 + ; (^-2) 

\ = Zi + f S2,g(Zi + Z2 — 2Z345) + S2i3C2i3{Zq — Z-j) , for Z = 3 , 4 , 5 , (A. 3 ) 

where Z 345 = Z 3 + Z 4 + Z 5 . However, these results do not take into account the fact that 

Ae = Ay = 0, which yields two relations among the Zj. These relations were given in Eqs. (26) 
and (27) and are repeated below for the convenience of the reader. Recall that we employ a 
convention where 0 < /5 < | 7 r. Then, Z 2 and Z 345 are dependent quantities for /3 7 ^ 0, ^tt, ^tt, 

Z 2 = Zi + 2(Z6 + Zy) cot 2/3 , (-^-4) 

Z 345 = Zi + 2Z6 cot 2/3 — (Ze — Zy) tan 2/3 , (A.5) 


An alternative form of Eq. (A.5) is obtained by combining the results of Eqs. (A.4) and (A.5), 
which yields 

Z 345 = Z 2 — 2 Zy cot 2/3 — (Ze — Zy) tan 2/3 . (A. 6 ) 

Taking the average of Eqs. (A.5) and (A. 6 ) provides one more useful relation that can be used 
as the second condition for the softly-broken Z 2 symmetry along with Eq. (A.4), 

■^345 = \{.Z\ -\- Z 2 ) -|- 2(Z6 — Zy) cot 4/3. (A.7) 

Using Eqs. (A.4) and (A.7) it follow that if /3 = 0, then Ze = Zy = 0 ; if /3 = | 7 r, | 7 r then 

Z 345 = |(Zi -|- Z 2 ); and if /3 = ^TT then Zi = Z 2 and Ze = Zy. 

Consequently, the expressions for the A* in terms of the Zj can be written in numerous 

equivalent ways depending on the choice of independent quantities. Eor example, if /3 7 ^ ^tt, 
then eliminating Z 345 and either Zi or Z 2 yields 


Ai 


A2 


Zi — Zg tan 2/3 -|- I tan^ /3 tan 2/3(Z6 + Zy), if /3 7 ^ |7r , 

Z 2 + Zy tan 2/3 — I cot^ /3 tan 2/3(Z6 -|- Zy), if /3 7 ^ 0 , 

Zi — Zg tan 2/3 -|- I cot^ /3 tan 2/3(Z6 -|- Zy), if /3 7 ^ 0 , 

Z 2 -|- Zy tan 2/3 — I tan^ /3 tan 2/3(Zq + Zy), if /3 7 ^ , 


(A.8) 

(A.9) 


Aj = Zj-|-^(Zg — Zy) tan2/3 , fori = 3,4, 5. (A. 10) 

Note that the Z 2 -basis and the Higgs basis coincide if /3 = 0 (in which case •hi = Hi and 
'I >2 = H 2 ) and if /3 = ^TT (in which case $1 = H 2 and 'I >2 = Hi). The two alternative forms 
given in Eqs. (A.8) and (A.9) are a consequence of the symmetry of Eqs. (A.4)-(A.7) under 
the interchanges, Zi < —Z 2 , Ze < —Zy, (3 i —— /3. 
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The exclusion of /5 = in Eqs. (A.8)-(A.10) is an artifact of expressing these results in 
terms of both Zq and Z-j. Nevertheless, there is no discontinuity, since Zq = Zt Ski (3 = 

One way to avoid this inconvenience is to eliminate either Zq or Z^ in favor of ^ 345 . The end 
result is 



Zi(l 

— 1 tan^ 13) -b 

IZ345 

tan^ P - 

■ ^Zg tan/ 3(5 ■ 

- tan^ P ), 

lip 

7^ \tt 

Z2(l 

— ^ cot^ /3) + 

IZ345 

cot^ P — 

^Z7 cot / 3(5 - 

- cot^ P ), 

up 

7^ 0, 

Zi(l 

— 1 cot^ /3) + 

IZ345 

cot^ P -b 

|Z 6 Cot/ 3(5 - 

- cot^ P ), 

up 

7^ 0, 

Z2(l 

— 7 tan^ P) + 

^Z 345 

tan^ P -f 

- ^Z7 tan/ 3(5 - 

- tan^ P ), 

if P 

7^ Itt 

Zi + 

|(Zi — Z345) - 

T Zq cot 2/3, 

for / = 3 , 4 , 5 

and /3 7^ 0, 

Itt, 


Zi + 

i(Z2 — Z345) ■ 

- Z7 cot 2/3, 

for / = 3 , 4 , 5 

and /3 7^ 0, 

iyr. 



(A.11) 

(A.12) 

(A.13) 


Finally, one may choose to eliminate both Zq and Z 7 , using Eqs. (A.4) and (A.7). The end 
result is valid for /3 7 ^ | 7 r, ^tt, 

1 

Ai = i(Zi + Z 2 ) + -^(Zi + Z 2 - 2 Z 345 ) + ^(Zi - Z 2 ), (A.14) 

4C4/3 ^C2j3 

1 

A 2 = i(Zi + Z 2 ) + -^(Zi + Z 2 - 2 Z 345 ) - ;^(Zi - Z 2 ), (A.15) 

4c4/3 2c2/3 

^2 

Ai = Zj — ^ (Zi + Z 2 — 2 Z 345 ), for i = 3,4,5. (A. 16) 

2 C 4 ^ 

The conditions for stability of the scalar potential [Eq. ( 1 )] for Ae = A 7 = 0 were first given 
in [58], 

Ai > 0 , A 2 > 0 , A 3 > —\/ A 1 A 2 , A 3 + A 4 — I As I > —\/ A 1 A 2 • (A. 17) 

Using the results of this Appendix, one can rewrite the stability conditions in terms of the Zj. 
The resulting expressions are not especially illuminating, so we will not exhibit them explicitly. 

In addition, we note that (under the assumption of Ae = A 7 = 0) the Aj {i — 1,2,..., 5) can 
be reconstructed in principle as follows. Assume that cp-a has been deduced from precision 
measurements of the SM-like Higgs boson (assumed to be h), and (3 is determined via the 
properties of the heavier Higgs states. We also assume that all four Higgs masses (m^,, niHi tua 
and mH±) have been measured. Lastly, we assume that a small deviation in the signal strength 
for /j —>■ 77 can be attributed to the presence of a charged Higgs loop,^^ in which case we can 
extract a value for ghH+H-- With all this information in hand, we begin by using Eq. (47) [or 
equivalently, Eq. (42)] to obtain Zq. Next, we employ Eqs. (41) and (43) to obtain Zi and Z 5 , 
and Eqs. (33) and (34) for the squared-mass difference, to deduce Z/^ — Z 5 , which 

together with the previous determination yields a value for Z 4 . Close to the alignment limit, 

^^Eliminating both Zq and Z 7 is not particularly useful in the cases of /? = | 7 r, | 7 r, where Zi + Z 2 = 2 Z 345 
and in the case of /3 = ^tt, where Zi = Z 2 [cf. Eqs. (A.4) and (A.7)]. 

^^In absence of a clear deviation from the SM in the 77 signal, one would be forced to seek out some measurable 
triple Higgs coupling involving no more than a single SM-like Higgs boson to avoid a suppression of the term 
that is sensitive to Z 3 or Z 7 [cf. Eqs. (59)-(66)]. 
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we can use ghH+H- to extract Z 3 [cf. Eqs. (65) and (76)]. We now have enough information to 
evaluate .^ 345 . Finally, we can use Eqs. (A.9) and (A. 10 ) to obtain Z 2 and Z 7 . We now have 
all the Zi (for i = 1,2,...7), which can then be employed with the formulae provided in this 
Appendix to obtain the Aj (i = 1 , 2 ,..., 5). 


B Trilinear Higgs self-couplings in terms of physical Higgs 
masses 


It is convenient to re-express the trilinear Higgs self-couplings in terms of the physical Higgs 
masses. First, Eqs. (32) and (34) yield 

(Z 3 -V Zi^ — Z^)v^ = ‘2{m\ — m^) -|- Z^v^ + 2 ZqV^ cot 2(3. (B.l) 

Using this result along with Eqs. (41)-(43) and (49), we end up with 

[(Z 3 Z 4 - Z^)sp_a + = [m\ + 2{m\ - rh^)]s/ 3 _« + 2cot2l3{m\ - m^)c/ 3 _„ , (B.2) 

[(Z 3 -F Z 4 - Z5)ci3-a - Zys^-JIv"^ = [m% + 2{m\ - m^)]c; 3 _a -2cot2/3{m% - .(B.3) 

Noting that Eqs. (33) and (34) yield m\ — rnjj± = ^{Z 4 — Z^)v^^ the above results immediately 
yield 

[Z^sn^oc + ^ 7 C/ 3 -a] = [ml + 2{m]j± - m^)] S/ 3 _„ -F 2 cot 2^(m^ - rn‘^)cfi_ec , (B.4) 

[Z^cp_a - = [m]j -h 2(m|± - m^)]c, 3 _« - 2 cot 2l3{m% - m‘^)s(s-a ■ (B.5) 

Thus, from Eqs. (59)-(66), we obtain 

9hAA = + 2(m^ - rn^)]s/ 3 _„ -h 2cot2^(m^ - m^)c, 9 _„| , (B. 6 ) 

9haa = -rh^)]c/ 3 _„ -2cot2P{m% -m^)s, 3 _„| , (B.7) 

9hHH = - 2m% -ml + 2(3m^ - 2m% - ml){s/^-a cot 2/3 - C; 3 _a)c^_„ I (B. 8 ) 


9Hhh = 2^«ft + 2(3m^ - m% - 2ml){ 8 / 3 -^ cot 2/3 “ C; 3 _a)c^_„| (B.9) 

9hhh = + 2(m^ - m^){Cfi_aCOt2/3 + S;3_„)c5_„| , (B.IO) 

Qhhh = - 2(m| - rh^)(s/ 3 _„cot 2 ,d - C/ 3 _„)sJ_„| , (B.ll) 

9hH+H- = -“I + 2(m^± - m^)]s/3-a + 2cot 2/5(m^ - m^)c; 3 _a| , (B.12) 

9hh+h- = “M ['^H + 2{m]j± - m^)]c; 3 _„ - 2cot2/5(m| - I. (B.13) 
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